Theoretical analysis of the experiments on the double-spin-chain compound — KCuCls 
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We have analyzed the experimental susceptibility data of KCuCla and found that the data are 
well-explained by the double-spin-chain models with strong antiferromagnetic dimerization. Large 
quantum Monte Carlo calculations were performed for the first time in the spin systems with frus- 
tration. This was made possible by removing the negative-sign problem with the use of the dimer 
basis that has the spin-reversal symmetry. The numerical data agree with the experimental data 
within 1% relative errors in the whole temperature region. We also present a theoretical estimate 
for the dispersion relation and compare it with the recent neutron-scattering experiment. Finally, 
the magnitude of each interaction bond is predicted. 
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The low-dimensional quantum systems with the exci- 
tation energy gap are now attracting much interest both 
experimentally and theoreticallytl Possibility of the high- 
Tc superconductivity upon doping carriers to the gapped 
insulator lies as the background. The spin-ladder model 
and its realization in a real compound SrCu203 |ffiay 
be the most well-known substance for this scenarioE Of 
course, various other new compounds are also synthesized 
which can be explained by the low-dimensional quantum 
spin Hamiltonian.B~@ These experimental achievements 
now offer many data left for theoretical investigations. 

Tanaka et alu have measured the magnetic suscepti- 
bility of the KCuCla, and proposed that it can be ex- 
plained by the double-spin-chain Hamiltonian. The sus- 
ceptibility data show the spin-gap behavior. They es- 
timated the amount of the excitation gap A by fitting 
the low-temperature data with its theoretical expression, 
X{T) (X T"i/2exp(-A/T), given by Troyer et alE The 
estimated gap, A/fce 35K, is consistent with their re- 
cent measurement of the magnetization ^process at low 
temperature, which gives A//cb ~ 31.1K.E3 

For the thorough understandings of the system, we are 
necessary to determine the strength of each interaction 
bond of the model Hamiltonian which expresses the sub- 
ject material. For frustrated systems, one usually cal- 
culates all the eigenvalues of a given model Hamiltonian 
on a finite lattice by the numerical diagonalization, and 
then compares the obtained physical quantities at finite 
temperatures with the experimental data. Since the nu- 
merical diagonalization technique is restricted to the very 
small sizes with sixteen S = 1/2-spins or the less, the re- 
sults suffer a severe finite size effect. The quantum Monte 
Carlo (QMC) method can handle much larger systems, 
but has not been applied to the frustrated systems be- 
cause of the negative-sign problem. Recently, one of us 
(T. N.)|ii|as solved this problem in the double-spin-chain 
system.O It makes possible to treat large system sizes 
even at very low temperatures, and thus we are able to 
compare directly the numerical results to the experimen- 
tal data without suffering the size effect. 



In this Rapid Communication, we calculate the suscep- 
tibility of the double-spin-chain model with 162 spins, 
and aim at the determination of the interaction bonds 
that explain KCuCls. For this purpose, it is necessary 
to calculate the dispersion relation of the excited state 
to compare with the neutron-scattering experimental re- 
sults. The susceptibility alone is not sufhcient to dis- 
cuss the details of the system, since it is the integrated 
quantity. We use our analytic expression for the disper- 
sion relation obtained recently for general double-spin- 
chain models £3 Our analysis partly disagrees with the 
experiments, which may be understood by possible two- 
dimensional couplings. 

We consider the following spin Hamiltonian under the 
open boundary conditions: 
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Here, A'" is the linear size of the system, and |cr| = \t\ = 
1/2. Figure ^ shows the shape of the lattice. 
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FIG. 1. Shape of the general double spin-chain model. 

Three cases of the above Hamiltonian, (i) J\ = J3, (ii) 
J3 = 0, and (iii) Ji = 0, are especially investigated. In 
each case, J2 is set variable ranging from J2 = —8 to 
J2 = 8. The first one, Ji = J3, is what we call here the 
'zig-zag' or commonly called the railroad-trestle model. 
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This model has frustration for J2 > 0, and thus has not 
been analyzed on its thermodynamic properties yet. The 
second one is the ordinary two-leg 'ladder' model, and 
the third one is the 'bond-alternation' model. By chang- 
ing the sign of the J2-bond, the Hamiltonian can express 
both the dimer-gap system and the Haldane-gap system. 
Difference in the origin of the gap affects the structure 
of the excited states and thus the finite temperature be- 
havior of various physical quantities. We also comment 
on this point at the end of this Rapid Communication. 

Before the demonstration of the numerical evidences, 
we briefly summarize the simulational technique. We 
have done the ordinary world-line QMC simulations, but 
the choice of the representation basis is considered dif- 
ferent from the conventional one. Two spins, cr„ and 
T„, are coupled and considered as a unit of update. This 
dimer unit takes four states associated with the eigen- 
values of each spin, |(t^,t^). We rearrange these four 
states so that they have the spin-reversal symmetry: 

«i-(IT,T> + U,i»/V2, (2) 

«2 = (|T,T)-Ua))/x/2, (3) 

^'3 = (IT,i> + U,T»/V2, (4) 
^'4-(|T,i>-U,T»/V2. (5) 

Here, t and .[ denote the eigenstates. Now, two Ji- 
bonds and one Js-bond just become a single effective 
bond connecting the neighboring dimer units, and a J2- 
bond only contributes to the inner energy of a dimer unit. 
We can remove the negative-sign problem by using this 
representation basis. The Trotter number m are chosen 
so that /3{2Ji + J3)/m varies from 0.35 to 0.15 at each 
temperature, and the extrapolations to the infinite m are 
done. Typical number of the Monte Carlo steps is 500 
000 divided into ten parts to see the statistical deviation. 
The first 50 000 steps are discarded for the equilibrium. 
The autocorrelation time of the susceptibility is less than 
an order of unity. 

We have performed simulations for several values of J2 
in each choice of the Ji- and the J3- bonds. First, we 
have set Ji = J3 and done simulations thoroughly by 
changing a value of J2 from —8 to 8 in a step of 2. A 
rough estimate for the value of J2 is made at this stage. 
Then, we have proceeded to the other cases, Ji = and 
J3 = 0. 

Figure ||(a) shows the QMC results of the susceptibility 
compared with the experimental data, and (b) is its low- 
temperature plot. The experimental data are the mini- 
mum susceptibility when the magnetic field is applied in 
the cleavage plane. The absolute value of the susceptibil- 
ity depends upon the direction of the field, which can be. 
understood by considering the anisotropy of the g-value.cl 
Here, the 5-value is determined by the comparison at high 
temperatures, T > lOOK. Error bars of the QMC data are 
almost negligible. Magnitude of each interaction bond is 
determined in order that the low-temperature data agree 
quantitatively with the experiment. 
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FIG. 2. (a) The QMC results of the susceptibility (sym- 
bols) compared with the experimental data (line); (b) its 
low-temperature plot. The number of the spins is 162 and 
the open boundary conditions are used. The strength of each 
interaction bond is; (A) Zig-zag: Ji — = 8.35K, J2 — 6J1, 
g = 2.05, (o) Ladder: Ji = 12.3K, J3 = 0, J2 = 4Ji, 3 = 2.05, 
and (x) Bond-alternation: J3 = 12. 9K, Ji = 0, J2 = 4J3, 
g = 2.03, respectively. 

Triangles are the best fit to the experiment so far in 
the case of Ji = J3, and the values of the interaction 
bonds are: Ji — J3 = 8.35K, J2 = 6J1, i.e., all the 
interaction bonds are antiferromagnetic. The g-value 
g — 2.05. Circles are those of the ordinary ladder case, 
and the values are: Ji = 12. 3K, J3 = 0, J2 = 4Ji, and 
g = 2.05. Crosses are the case of the bond-alternation 
model: Ji = 0, J3 = 12. 9K, J2 = 4Ji, and g = 2.03. All 
three data quantitatively agree with the experiment at 
all the temperatures. Our estimate for the g-vahie is alsa 
consistent with the ESR- measurement giving g = 2.05.1121 
Relative errors from the experiment are within 1% in 
every case, and are smallest in the 'ladder' case. How- 
ever, we are not sure the small deviations among three 
cases are relevant or not, because we have used a very- 
simplified model spin-Hamiltonian for a real compound. 
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Systematic errors from an adoption of this model may be 
the most important one. Therefore, we cannot determine 
the strength of each interaction bond from the fitting 
of the susceptibihty data alone. Only a common fea- 
ture known from this plot is that the J2-bond is strongly 
dimerized antiferromagnetically, and thus the origin of 
the gap is the dimer gap. 

We try to determine the interaction bonds by the dis- 
persion relation of the excitation energy now. This value 
is quite sensible to the details of the model. Recently, 
Kato et al.tB applied the neutron-scattering analysis on 
this compound, and presented the dispersion relation. 
We compare our analytic dispersion relation with this 
experiment. 

In our previous paper ,0 we have deduced an expression 
giving the dispersion of the general double-spin-chain sys- 
tems as a function of the interaction bonds. There, we 
employed a single domain-YaU-ivariation after the non- 
local unitary transformationti'Ej is applied. This approx- 
imation is quite excellent near k — t:, checked with the 
numerical diagonalization results of finite systems, and 
with the exactly-known results. On the other hand, it 
becomes overestimated as the wave number k approaches 
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FIG. 3. The dispersion relation of the excited states of the 
system. Zig-zag: Ji — J3 = 8.35K and J2 = 6J1, Ladder: 
Ji = 12. 3K, J3 = 0, and J2 = 4Ji, and Bond-alternation: 
J3 = 12. 9K, Ji = 0, and J2 = 4J3, respectively. The neu- 
tron-scattering experimental data are plotted by circles for 
the constant wave number perpendicular to the chain k± = 0, 
and by solid squares for k± = tt. 

By using the values of the interaction bonds obtained 
from the susceptibility fitting, we give the dispersion of 
the excitation energy in Fig. ^. The first excitation is 
at A; = TT in the 'zig-zag' and the 'ladder' models, and 
is at A: = in the 'bond-alternation' model. Within 
our analysis, the first excitation is always at k = tt for 
Ji > min( J2, J3)/2, and is otherwise at k — 0, where 
min(J2, J3) stands for the minimum value of J2 and J3. 

The dispersion of the 'ladder' case is consistent with 
the experimental dispersion parallel to the double chain 



when the constant wave number perpendicular to the 
chain k± is equal to zero. Our estimate of the gap is 
A/fce ^ 38K, which is a little larger than the results of 
the magnetization process experiment.E3 So the 'ladder' 
model may be the most favorable candidate for explain- 
ing the KCuCla at the present stage. However, it should 
be noticed that the experimental dispersion is dependent 
upon k±. If a double-chain is isolated from each other, 
there should not be the k± dependences. Therefore, it 
might be an evidence of two-dimensional interaction cou- 
plings. This poiat is left for the future study. 

Tanaka et als fitted their experimental data of the 
susceptibility to the theoretical expression -X{T) = 
{l/2^/naT) exp(— A/T), given by Troyer et alB suppos- 
ing the form — A+a\k—Tr\'^ for the magnon dispersion. 
This expression is valid when T <C A, Z?, where D is the 
band width of the magnon dispersion. If we take the 
effect of the band width into account, we obtain 



(6) 



where 7(2;, p) is the incomplete gamma function defined 

by 



(7) 
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Since 13 <C A from Fig. |l x(T) behaves as 



xiT) 




{T<t:D<t:A), 

(D < T< A). 

(8) 



Thus, the expression of Troyer overestimates the suscep- 
tibility. In the fitting by Tanaka, the theoretical curve 
of Troyer severely deviates from the experimental curve 
for T > 20K, which is well-explained by D = 2.2 meV = 
25 K of the 'ladder' model. 

In the last part of this Rapid Communication, we men- 
tion the difference of the susceptibility in the dimer-gap 
system and that in the Haldane-gap system. Only the 
'zig-zag' case, Ji = J3, is demonstrated as an example 
in Fig. ^. Product of the susceptibility and the tem- 
perature, xT, is plotted against the temperature, since 
this value is independent from the scale of the interaction 
bond. We use the logarithmic scale for the temperature 
axis, so that the rescaling of the temperature by chang- 
ing the values of interaction bonds just causes a parallel 
shift along the temperature axis. Circles are the data 
of the dimer-gap system: = J3 = 8.35K, J2 = 6J1, 
g — 2.05, and crosses are those of the Haldane-gap sys- 
tem: Ji ^ J3 = 45K, J2 = -4Ji, g = 2.00. These in- 
teractions are determined by fitting the low-temperature 
data, while the g-vahie is estimated at high tempera- 
tures. As seen from this figure, function form of xT in 
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the Haldane-gap system is obviously different from the 
experiment, so we cannot fit them by any parallel shift. 
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FIG. 4. Product of tlie susceptibility and tfie temperature, 
XT, is plotted against the temperature. The QMC results of 
the dimer-gap system (o) and the Haldane-gap system (x) are 
compared with the experimental data (line). The strength of 
each interaction bond is, (o) : Ji = J3 = 8.35K, J2 — 6J1, 
g = 2.05, and (x): Ji = J3 = 45K, J2 = -4Ji, g = 2.00, 
respectively. 

Let us consider the dimer limit, Ji = J3 = 0, J2 > 0, 
where the ground state is an array of independent sin- 
glet dimers and thus is non-degenerate. When the in- 
teractions between singlet dimers (Ji and J3) are intro- 
duced, the ground state is somewhat modified, but still 
has the nature of the singlet dimer. On the other hand, 
the ground state in the Haldane limit, Ji = J3 = 0, 
J2 < 0, is an array of independent triplet dimers and 
is highly degenerate. When Ji and J3 are switched on, 
the degeneracy is lifted up resulting in the unique ground 
state. The density of states of low- lying excited states are 
larger than the singlet dimer case, reflecting this high de- 
generacy of the ground state at Ji = J3 — 0. Therefore, 
the susceptibility peak in the Haldane system becomes 
lower and broader than in the dimer system. We may 
determine the origin of the gap from the full width at 
half maximum of the susceptibility, or the function form 
of xT. 

In summary, we have analyzed the susceptibility and 
the neutron-scattering measurements on the KCUCI3 
from the theoretical point of view. Large-scale quan- 
tum Monte Carlo simulations clarified this compound 
is well-explained by the double-spin-chain Hamiltonian 
with strong antiferromagnetic dimerization. Since the 
susceptibility alone is not sufficient to determine all the 
interaction-bond strengths, we also calculated the dis- 
persion relation of the excited state to compare with the 
recent neutron-scattering experiment. As far as we re- 
strict our theoretical analysis to the isolated double-spin- 
chain Hamiltonian, the best one to explain both measure- 
ments at present is the ladder model, which is defined by 
Ji = 12. 3K, J2 = 4Ji, J3 = 0, and g = 2.05. This g- 



value agrees with the ESR- measurement However, the 
experimental dispersion relation has strong dependence 
on the constant wave number perpendicular to the chain, 
fcx, which cannot be explained by our theoretical anal- 
ysis on the isolated double-spin-chain model. Thus, the 
system may have two-dimensional interactions that are 
not relevant at low temperatures. 

The quantum Monte Carlo simulation has become the 
most realistic method to investigate the various quasi- 
one-dimensional compounds that can be expressed by the 
general double-spin-chain Hamiltonian even if the system 
has frustration. Present calculations can be extended 
to analyze other experiments easily. For example, we 
comment on that of Cu2(l, 4-Diazacycloheptane)2Cl4 re- 
cently done by Hammar et aln Their susceptibility data 
are consistent with our calculation down to the lowest 
temperature that could not be obtained by the numer- 
ical diagonalization. The choice of the interactions is 
Ji = J3 = 2.26K, J2 = 6J1, and g = 2.04. 
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